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Adiabatic fusion potentials, extended well into compound nucleus region, are calculated for a 
number of reactions within the static Hartree-Fock method with the Skyrme SkM* interaction. 
The calculated fusion barriers agree quite well with the data, in spite of considerable errors in 
reaction Q values. This suggests some error cancellation, possibly with the relative kinetic energy 
term. Calculated barrier heights are consistent with the idea of fusion hindrance in tip collisions. 
Adiabatic potentials differ considerably from the results of the frozen density approximation. We 
briefly discuss relative positions and heights of the fusion and fission potentials and speculate on 
their possible connection with the fusion hindrance for large ZrZp. 
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INTRODUCTION 



Entrance channel potential is one of the principal ingredients in the description of heavy ion collisions. It controls 
not only the capture probability, but is also related to a substantial hindrance of the compound nucleus (CN) forma- 
tion, observed experimentally as a large probability of quasifission in reactions between targets and projectiles with 
' ZtZp >1800. In this way it is crucial for prospects of the synthesis of the heaviest elements. 
, Up to now, calculations of the entrance channel potentials were performed mostly within some phenomenological 
fSl ' models. The mean field Hartree-Fock (HF) methods were applied only with the assumption of frozen densities (see 
. Q Q for recent examples), which involves approximations spoiling selfconsistency. The same assumption underlies 
' the nuclear matter approach, see e.g. 0. It was sometimes argued that such instantaneous, diabatic fusion barrier 
should be expected in actual heavy ion collisions due to the short time scale involved in passing over the barrier. 
Nevertheless, too high barriers and unrealistic, strongly repulsive potentials obtained in the frozen density regime 
^ ^] point to a limited relevance of this method. We think that the adiabatic potential is a necessary ingredient in 
^ . the selfconsistent study of nuclear fusion in much the same way as the static barrier is a necessary first step in a study 
^ ' of nuclear fission. 

(— I , In this work, we have performed HF calculations for a number of target-projectile combinations and found the 
adiabatic nucleus- nucleus potential and (outer) fusion barrier. As might be expected, the adjustment of the nuclear 
matter to the mutual interaction of target and projectile gives markedly different results from those of the frozen 
density approach. 

5— i ■ Within the mean field, for a deformed target or projectile, there are many fusion barriers, depending on the relative 
orientation of fragments. We have calculated potentials for tip and side collisions, i.e. for the configuration with 
the symmetry axis of the deformed fragment parallel and perpendicular to the relative distance vector. Effects of 
the fragment orientation and deformation on fusion were proposed in Quite recently, experimental evidence was 
presented showing that at least part of the CN formation hindrance is related to the dominance of quasifission 

in tip collisions. 

We have used mostly the Skyrme SkM* interaction SI, originally invented to properly fit the fission barriers. We 
are aware of previous HF calculations of the potential energy in a dinuclear regime, like e.g. (see also recent work 
|lO|). however, their objective was nuclear scission, where deformations of separating fragments much differ from those 
expected in the fusion entrance channel. 

The aim of this paper is twofold. First, we provide some systematics of the selfconsistent fusion barriers and 
check against the experimental data and the recent calculations within the frozen density regime Q with the same 
interaction. In our study, we include some systems with well measured fusion cross sections as well as reactions on 
^"^Pb target of the typ e studied in GSI Darmstadt 0, 0| and ^*Ca-induced reactions used in recent experiments 
at JINR in Dubna [Hill El, both reported as leading to the synthesis of the heaviest elements. Second, we look 
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for a relative position (deformations) and height of the fusion channels with respect to the fission valley to figure out 
quantum transitions implied by the fusion process. Our original motive was to find whether the adiabatic energy 
surface offers some clue to the fusion hindrance. 

Some preliminary results of this study have been published in . Differences between the present barriers and 
those in [l&J come from more extensive HF iterations and the improved numeric treatment of kinetic energy. 

II. METHOD OF CALCULATIONS 

We calculate static potentials, which means that we neglect currents and other time-odd quantities in the Skyrme 
energy functional and consider time-reversal invariant HF states. This seems reasonable at energies close to the 
Coulomb barrier. The potential between nuclei 1 and 2 in the central collision is then 

V{R) = E{R) + Si + ^2, (1) 

where E{R) is the (negative) HF energy of a dinuclear complex at the distance R and Bi are the (positive) binding 
energies of target and projectile. In order to have a numerically consistent treatment, Bi, i = 1,2, and E{R) are 
calculated with the same HF code. In the present calculation pairing is neglected. 

The HF equations have been solved on a spatial mesh of a size proper to the colliding system. Our code assumes 
two plane symmetries, i.e. it allows for the mass asymmetry along one direction. We consider tip and side collisions, 
with the angle between the symmetry axis of a deformed nucleus and the line connecting centers of two fragments 
equal to 0° and 90°. 

Initially, two sets of wave functions corresponding to two fragments are placed at a chosen distance being an 
integer multiple of the mesh spacing (0.69 or 0.77 fm). The HF equations are solved by the imaginary-time evolution. 
Wave functions are kept orthonormal and this enforces the Pauli principle. For fragments placed close enough, the 
necessary rearrangement of orbitals occurs already at the beginning of the HF procedure and avoids higher than 
normal densities. Final wave functions correspond to the local minima of the energy functional to which the initial 
configuration converged. For compact configurations of the CN, i.e. small R, there may be quite a few such minima, 
corresponding to different valleys in the nuclear energy (hyper-)surface. Thus, in contrast to the fusion barrier, the 
entrance channel (fusion) potential V{R) is not unique. The distance R between two fragments is calculated as the 
distance between the mass centers of two half-spaces containing Ai and A2 nucleons. 

Since HF states are superpositions of various total momentum eigenstates, one has to include the center of mass 
(cm) kinetic energy correction. The expectation value of the cm kinetic energy for an A^-fermion Slater determinant 
reads 

N N 

Ecm{N) = ^^(E(^ I p2 I fc) - E I (fc I p I I) n, (2) 

k=l k^l 

with k, I labelling occupied single particle states. Only exchange two-body terms are present in Eq.(|21) as the diagonal 
matrix elements of momentum vanish in the time-reversal invariant states. In most of the HF calculations two-body 
terms in Ecm are neglected (for a competent review of the cm corrections in HF see (3)- With the Skyrme SkM* 
force, Ecm is taken as the average kinetic energy, (t) — J2k-occ(^ I ^ I fc)M- 

Independently of the approximations made, for a Slater determinant describing two separated fragments with A^i 
and iV2 nucleon wave functions, A^i -|- N2 = N, with vanishing mutual momentum overlaps, {ki | p | ^2) = 0, one 
obtains the relation 

EUN.) + EUN.) = EUN) + ^^^^-W) + A^i^c^(iV.) _ 

This implies that in any HF calculation which includes Ecm, the nucleus- nucleus potential V{R) of Eq.(^ will tend 
to the value of the second term on the right hand side of Eq.(|3J, when R tends to infinity. This term is just the 
expectation value of the kinetic energy of the relative motion of the two fragments. In order to preserve the usual 
meaning of the Coulomb barrier one has to subtract this asymptotic term in Eq. Q . Specifically, with the SkM* force 
and for fragments with Ai and A2 nucleons, Ai + A2 = A, one has to subtract (A2(ii) + Ai{t2))/A. 

This subtraction has been included in the frozen density calculations However, it is surely incorrect when 

two fragments become one, i.e., for compact configurations of the system. As we do not know how kinetic energy 
of the relative motion of the fragments transforms into potential energy of the combined system, and still want to 
compare with the existing frozen density results, we also calculate V{R) normalized to zero at infinity. One has to 
keep in mind that such potential underestimates the true V{R): the smaller the distance R, the larger the offset which 
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reaches —14 to —18 MeV at the CN ground state. The product ZtZp controls the overlap of target and projectile 
at the fusion barrier: the larger Z^Zp, the larger the overlap. Hence, one can expect that Ecm correction, valid for 
separated fragments, produces too low fusion barriers for fragments with large ZtZp. 

III. RESULTS AND DISCUSSION 

Examples of entrance channel potentials are shown in Figs. 1,2. In Fig. 1 we show the tip- and side collision barrier 
for the ^^^Pu-f^^Ca reaction. Two results are shown, differing by the number of performed HF iterations. Adiabatic 
(lower) barriers were obtained after three long series of iterations, required in view of an inefficient convergence of the 
pure HF without pairing. Only one series of iterations produced the other (higher) plotted barriers. The quantities 
obtained in this way were in fact published in |0| . It is seen that an incomplete convergence in the true HF mimics 
to some extent the frozen density results 0,0: for systems with large ZtZp there is a potential pocket inside the 
fusion barrier and V{R) rises steeply with decreasing R. In the frozen density method this results from a higher 
than normal density in the overlap region. In the HF, it is the initial orthogonalization which disturbs density in the 
overlap region and rises energy. A subsequent quenching is necessary to lose this excitation. 

In quantitative terms, even non-fuUy adiabatic selfconsistent potentials are lower and rise much less steeply with 
decreasing R than the frozen density results. This rise vanishes in fully adiabatic potentials, as seen in Figs. 1,2. These 
potentials decrease with decreasing _R, except for the superheavy systems, where the Coulomb barrier becomes a local 
plateau at the level of, or even below, the ground state of CN. In addition to the ^^^Pb+^^Zn reaction, potentials with 
a plateau in the place of the barrier are obtained also for the ^"^Pb+^^Ge and ^^^Sn+^^^Sn systems. A comparison of 
fusion barrier heights shows that for reactions on ^''^Pb and actinide targets leading to very heavy nuclei the adiabatic 
barriers are by 10-20 MeV lower than the frozen density ones 0] (see also Table 1). Taken at the face value, the 
adiabatic results suggest that, after passing the fusion barrier, systems with sufficiently low ZpZt experience a steady 
drivig force towards CN formation, while with increasing ZtZp this force weakens and even changes to repulsion. 

It may be noticed in Fig. 1 for the reaction ^''^Pu-|-'**Ca that the shorter and higher side collision potential offers 
larger driving force towards CN than the longer but lower tip collision one. The difference in the barrier position for 
tip- and side collision for ^'*^Pu-|-''*Ca system, Rup — Rside, is about 2 fm. Similar situation occurs for the reaction 
^^"Pd+^^^Pd (Fig. 2), for which the tip-tip collision potential shows a plateau with fluctuations. On the other hand, 
the higher side-side collision barrier gives chance to a quite steep (20 MeV) slide towards CN pocket (see below). The 
potential for the ^"^Pb-f ^"Zn reaction, on the other hand, does not provide any force in the CN direction. This could 
suggest that the entrance channel potentials give a better chance of fusion in the side collisions with prolate deformed 
actinide targets as compared to reactions on ^°®Pb target. 

When comparing various results for barriers one has to remember that with rising ZtZp 1) the frozen density 
approach becoms poorer, 2) the correction for cm energy becomes questionable and may produce too low barriers, 3) 
the adiabaticity effect on the lowering of the barrier increases. 

Calculated fusion barriers Bcai, taken as the locally highest value of V{R), rounded to 0.5 MeV, are compared in 
Table 1 to the Bass interaction barriers ^TS\ and to the recently given threshold barriers Bthre ■ The interaction 
barriers are chosen here instead of the fusion barriers since they much better approximate the experimental 
barriers for heavy systems. The threshold barriers [l9l | are derived from sufficiently detailed fusion data and should 
correspond to the calculated adiabatic barriers. The values of Bthre for the heaviest systems are based on the capture 
data [2^. For deformed targets, both the calculated tip and side (in parentheses) collision barriers are given. In 
order to appreciate excitation of a compound system at the fusion barrier we also give in Table 1 the experimental 
(or extrapolated) reaction Q values. 

Out of nine threshold barriers, two are overestimated and six underestimated in the calculations by up to 4 MeV 
(it seems that Bthre for the ^^^U-|-^^0 reaction corresponds to tip collision |^). One, for the ^•^*U-|-^^Ca reaction, is 
6 MeV below the calculated side collision barrier and 11 MeV above the tip collision one. Although for two reactions 
leading to ^^^Th CN the barriers Bthre are not known exactly, their rough estimate from data [ij is possible. For the 
reaction i^2^_j_32g ^j^g estimate agrees well with our Bcai, while for ^^^Sm-|-^°Ni it seems close to the calculated tip 
collision barrier, although the lack of the low energy data causes large uncertainty. It should be stressed that the lack 
of observed evaporation residues (ERs) at Ec„i=175 and 182 MeV in the experiment served as an argument for 
the fusion hindrance in tip collisions. 

The fusion barrier calculated for the recently measured reaction ^^^Sn-|-^'*Ni is higher than Bthre suggested by the 
data |(25|, but lower than the Bass interaction barrier assumed in the coupled channels analysis there. The ^^°Pd 
nucleus, previously considered vibrational (see e.g. H^), is now well established by the Coulomb excitation studies 
|26| as a rotor with deformation (32 ~ 0.25. We have obtained a shallow prolate-deformed HF minimum at Q =7.4 b, 
consistent with this, and used it to calculate tip-tip and side-side collision barriers given in Fig. 2 and Table 1. The 
^^°Pd-|-^^°Pd— >^^°U reaction is a very well known case of the severe fusion hindrance |23, at least by four orders of 
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magnitude. The lowest measured point at Ecm =228 MeV is already above the calculated side collision barrier. 

For most reactions leading to superheavy nuclei there are too few data to extract threshold barriers, even from 
capture data. However, the data on evaporation residue formation give some upper bound on the height of the fusion 
barrier, as their production at Ecm below the lowest barrier is improbable. For ^°^Pb target, the ^^^110 ERs were 
detected following reaction with "^"^Ni at Ecm = 235.3 - 238.2 MeV JJtiM^ while in reaction with ^"Zn two observed 
^^^112 ERs were produced at Ecm. =257.2 and 259.1 MeV In reactions with ''^Ca projectiles on actinide targets, 
two ERs were observed in the reaction on ^ss^j _ 190.5 — 193.9 MeV three events in the reaction on 

244pu target at Ecm = 194.5 - 202 MeV [Hll^l and one event for '^'^^Cm target at Ecm = 199.7 - 205.1 MeV 
All these bombardment energies in the cm. system are higher than our calculated barriers: by at least 6 and 13 MeV 
for ^^Ni and ™Zn projectiles on ^'^^Pb target, and by 2-3 MeV higher than the calculated side collision barriers for 
actinide targets. 

Overall reasonable agreement of the calculated fusion barriers with data is quite astonishing in view of the deficiency 
shown by the SkM* force in reproducing the experimental binding energies. As follows from our calculations and the 
recently calculated mass table |3ll | there is a large underestimation of binding for heavy and superheavy nuclei with 
the SkM* force, e.g. about 9 MeV for ^^^W and ^'^'^Pu, but aheady about 14 MeV for ^sgnq, and expected 18 MeV 
for ^^^114, according to usually fair estimates of the phenomenological Thomas-Fermi model Moreover, there 

are some errors specific to SkM* force, like too large bindings of ''^Ca (by 4 MeV) and ^'^^Sn (by 7.5 MeV). 

It is reasonable to ask about a correlation between errors in fusion barriers and those in reaction Q values. It turns 
out that the Q values are overestimated by up to 5 MeV for lighter compound systems, i.e. similar in magnitude 
but mostly opposite to the errors in barriers. However, for heavier compound systems the errors in barriers are much 
smaller than those in Q values, e.g. Q values are too large by about 11 MeV for ^^Zr-l-^^Zr, about 17 MeV for 
"'^Ca-l-^^^Pb, expected about 13 MeV too large for ^^Ca-|-^'*'*Pu. All this points to some error cancellation, which, as 
discussed in the previous section, might be related to the overestimated relative kinetic energy correction. 

An intriguing question is whether the calculated potentials V{R) may be related quantitatively to the fusion 
hindrance seen in experiment. In P|, a sufficiently deep pocket in the frozen density potential V{R) was proposed as a 
necessary precondition for fusion. The quantification of this assertion is difficult, as the depth and the very existence 
of the pocket relics on the unrealistic frozen density intrinsic barrier, which exists e.g . already for ^"Zr-l-^^Zr and 
208pj^_l_48Q^ reactions, for which no fusion hindrance is experimentally observed |20ll2]| . 

Reasoning intuitively, a relative situation and heights of fusion and fission valleys should be relevant for prospects of 
fusion: to be captured in the CN configuration, a system must pass inside the CN fission barrier. In order to gain some 
orientation in this respect, a more detailed calculations have been performed for a few nuclei. The results for ^^^114 
^244p^_l_48Qg^ reaction) are shown in Fig. 3. When showing fission and fusion barriers in one figure the kinetic energy 
of the fragment relative motion has not been subtracted from the fusion barrier. Instead of the distance R between 
two fragments the quadrupole moment Q2Q of the compound system is shown on the abscissa. The calculated mass- 
and axially symmetric fission barrier of 15 MeV is heavily overestimated as compared to the experimental estimate 
Bfis ~ 6.5 MeV [s^, consistent with the Strutinsky method calculations [s^. This comes partly from the lack of 
pairing and neglect of nonaxiality. Accounting for both these effects still leaves a few MeV discrepancy between the 
data and the barrier obtained with the SkM* force. Similarly, the second hump of the barrier will be lowered by 
pairing and mass-asymmetry. Therefore, Fig. 3 gives a general idea about relative position of important points of the 
nuclear energy landscape rather than exact barrier heights. 

The tip-collision potential lies below the second hump of the fission barrier and corresponds to a sizable mass 
asymetry, increasing with Q2o- In the right portion of the figure, we extend fission barrier to mass asymmetry by 
showing cuts through the potential energy surface at fixed Q2Q =100, 125 and 140 b vs. octupole moment Q30. It is 
seen that the mass-asymmetric minimum becomes energetically well separated from the symmetric fission valley with 
increasing Q2q- The nuclear states building fusion potential are very close to the states in the fission valley with alike 
quadrupole and octupole moments, e.g. at Q20 = 140 b and Q30 ~ 25 • 10^ fm^, at Q20 ~ 125 b and Q30 w 20 • 10"^ 
fm'^ and at Q20 ~ 100 b and Q30 ~ 12 • 10^ fm^. Slight differences in energy (1-2 MeV) come from differences in 
necking. One can say that the tip coUission potential arises from traversing potential energy surface, in this case not 
quite along the conditional minima at fixed Q201 as these lie slightly lower. 

This overall picture suggests that, after surpassing the fusion barrier, the system has a good chance to fall inside 
the outer fission barrier. If some energy is lost to excitation on the way, it may be stopped at the higher inner barrier. 
One may expect that fusion in the tip collisions is hindered in this case. 

The side collision potential lies above both humps of the fission barrier and, after surpassing the fusion barrier, the 
system can roll down all the way to the CN ground state, even if some energy is transferred to excitation. (Looking 
at Fig. 3 it is good to remember that the side collision barrier has a large axial asymmetry.) 

For ^^^No (^"^Pb-t-^'^Ca reaction, not shown) the mass-asymmetric entrance channel potential leads inside the 
higher, inner hump of the mass-symmetric fission barrier. The calculated height of the latter, 15 MeV, again exceeds 
the experimental estimate consistent with the value w8 MeV obtained in the Strutinsky-like studies, but is similar to 
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the barrier (12.6 MeV) calculated for ^^''No with SkLy4 force [s^l- The lower, secondary hump of the fission barrier 
does not play any role. Therefore, one should not expect drastic fusion hindrance in this reaction. 

The above are examples of possible qualitative inferences one can make from knowing adiabatic energy surfaces. 
Surely, the concept that the relative situation of the fusion potential and fission barriers has relation to fusion hindrance 
requires further study. It is plausible that not only nuclear configurations, but also configuration changes which occur 
with some probability, are involved in the resulting loss of the probability of fusion. 

Conclusions of this work may be stated as follows: 

- Adiabatic fusion barriers calculated with SkM* force compare rather well with experiment, anyway much better 
than the frozen density results. 

- The barriers compare with experiment much better than the reaction Q values for heavy and superheavy systems. 
This may suggest that there is some error cancellation: subtraction of a too large relative kinetic energy makes up 
for too large reaction Q values. 

- The rise of the potential V{R) with decreasing R seen in the frozen density or not fully adiabatic calculations 
vanishes in adiabatic potentials, except for the heaviest CN where the fusion barrier top lies below the ground state. 
Thus, the relation of this rise to the fusion hindrance is not clear. 

- The comparison of our calculated barrier heights to data supports the idea that in reactions with large ZtZp and 
deformed fragments fusion is suppressed in tip collisions in spite of the substantially lower fusion barrier. 

- Judging from forces driving towards the CN minimum, the formation of superheavy CN may be easier in side collisions 
with deformed actinides used in JIHR Dubna, than in reactions with spherical targets used in GSI Darmstadt (without 
considering survivel probability). 

The latter two points suggest that the size of a necessary rearrangement of the nuclear matter drop is more vital 
for fusion than the barrier heights. One may wonder whether some measure of this rearrangement is not a proper 
variable to describe the fusion hindrance. One such variable could be Rfus — Rfis, where the position of the inner 
fission barrier is relevant. Obviously, these speculations call for a more detailed study, including multidimensional 
energy surfaces and pairing. 
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TABLE 1 - Calculated fusion barriers for tip (side) collisions in MeV vs. threshold and Bass interaction (not 
fusion) barriers 1I81I. Threshold barrier for ^°Zr+^°Zr is inferred from that for from J^L and for 

Pd+'^^S from |22|. R.eaction Q values (experimental/extrapolated or predicted (with asterisk) |2J|) are given 

in column 2. 



riv^f Pin 





-'-'cat 


^thre 


rfint 
^Bass 


40ca+40Ca 


14.3 


52 


5G.2±G.2 


52.5 




48.8 


86.5 




94.5 


90zr+40Ca 


57.3 


92 


92 7±0 6 


97.7 


96Zr+40Ca 


41.1 


85 




96.6 


902r+90Zr 


157.3 


175 


^^175. 85 


181.0 


"°Pd+32g 


35.4 


78f79 5"! 


80 4±0 2 


88.9 




199.7 


200f222 51 




228.3 




84.9 


121 5^1 24"! 




134.3 


i54Sm+*^0Ni 


147.6 


17Gfl85 5) 




191.5 


238U4.160 


38.3 


67f71) 


--71 


82.9 


i32Sn+64Ni 


111.1 


147 




155.7 


i32Sn+i32Sn 


260.8 * 


254.5 




257.4 


208pb+48Qj^ 


153.8 


170.5 


169±2 


176.1 


208pb+64Ni 


225.1 


229 




241.3 


208pb+™Zn 


244.9 * 


244 




256.3 


208pb+82Gg 


262.5 * 


260.5 




268.8 


238u+48Ca 


160.8 * 


171(188) 


182±2 


193.8 


244pu+48Ca 


163.0 * 


174(192) 




197.3 


248Cm+48Ca 


169.3 * 


183(196.5) 




201.0 


250cf+48Ca 


177.0 * 


187.5(201) 




205.1 



V(R) [MeV] 



200. 
190. 
180. 
170. 



x \ .a S-N 

x \ ^.-x ..,-B" 

13' % 



■•.a 



8. 10. 12. 14. 



R[fm: 



FIG. 1: Partially and fully adiabatic potentials for ^**Pu+**Ca tip (pluses, stars) and side (crosses, squares) collision obtained 
with Skyrme SkM* force. 
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FIG. 2: Adiabatic nucleus-nucleus potentials for the reactions: ^^Zr+^^Ca (shifted up by 40 MeV, open squares), ^^^Sn+^^Ni 
(full squares), ^°**Pb+^*'Ca (stars), ^°Zr+^°Zr (pluses), "°Pd-|-"°Pd tip-tip (open circles) and side-side colUsion (full circles) 

and ^"'^Pb+^Zn (crosses). 
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FIG. 3: Left panel: Mass-symmetric fission barrier (pluses) and fusion barriers for the "'^Ca-induced reaction (tip collision - 
crosses, side collision - stars) for the compound system ^^^114. Right panel: fission valley in the same system vs. octupole 
moment Q30 for Q20 = 100 b (stars), 125 b (crosses) and 140 b (pluses) . 



